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Spectroscopy with Multichannel Correlation Radiometers 



in 
o 
o 

(N 
u 

< 
O 



> ■ 
On : 

Tt- ■ 
^t- ■ 

o : 

in . 
o ■ 

Or 
i 1 

O ■ 

a: 



13 



A.I. Harris 

Department of Astronomy, University of Maryland, 
College Park, MD 2074%\ 

Correlation radiometers make true differential measurements in power with high accuracy and 
small systematic errors. This receiver architecture has been used in radio astronomy for mea- 
surements of continuum radiation for over 50 years; this article examines spectroscopy over broad 
band widths using correlation techniques. After general discussions of correlation and the choice 
of hybrid phase experimental results from tests with a simple laboratory multi-channel correlation 
radiometer are shown. Analysis of the effect of the input hybrid's phase shows that a 90° hybrid 
is likely to be the best general choice for radio astronomy, depending on its amplitude match and 
phase flatness with frequency. The laboratory results verify that the combination of the correlation 
architecture and an analog lag correlator is an excellent method for spectroscopy over very wide 
bandwidths. 



I. INTRODUCTION 



Modern analog lag correlators are capable of autocor- 
relation spectroscopy over wide bandwidths Q] . This ar- 
ticle examines an application of analog lag correlators, 
that of measuring cross-correlation functions for spec- 
troscopy with correlation detection techniques. Although 
variants of the correlation detection schemes are common 
for radio continuum radiometry, it seems that no one has 
yet adapted the architecture for spectroscopy (multichan- 
nel radiometry) for high-resolution spectroscopy with a 
single telescope. A cross-correlation spectrometer on a 
single-aperture instrument would have the same excellent 
stability as cross-correlation spectrometers in aperture 
synthesis arrays. In addition, as shown in Section llVl a 
multichannel correlation radiometer can share one ana- 
log cross-correlation backend spectrometer between two 
sky positions, providing a dual-beam system with about 
half the backend spectrometer cost and complexity of a 
receiver with dual total-power spectrometers. 



Correlation radiometers have made accurate radio and 
millimeter-wave continuum intensity measurements of 
the radio sky, the Cosmic Microwave Background, of 
rapidly changing scenes, and polarimetry (e.g. 0, y, S 
5, 6J). A number of authors have described specific ar- 
chitectures and examined the operation and sensitivity 
of correlation radiometers in absolute terms and their 
suppression of effects from the 1/ f noise common to am- 
plifiers 0, S H 1 E IH E El H . In the following, 
Section|n]contains a general discussion of correlation and 
Section lTTTl gives an analysis of the choice of hybrid phase, 
information that is not readily available elsewhere. Sec- 
tion llVl shows experimental results from tests with a sim- 
ple laboratory multichannel correlation radiometer that 
verify the theoretical expectations. 
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II. CORRELATION DETECTION 

Ryle introduced correlation detection to radio astron- 
omy with his invention of the phase-switching interferom- 
eter ^IJ. Ryle's interferometer squared the sum of the 
voltages from two antennas after modulating the phase 
of one antenna's signal. Maintaining phase sensitivity by 
multiplying voltages instead of detecting the total power 
alone made this a correlating instead of a phased array 
of antennas. Phase switching the signal from one an- 
tenna was the key element in the method's success, as it 
separated the desired cross-product of the two antenna's 
voltages from total power signals from the individual an- 
tennas. Correlation detection with rapid phase switching 
brought a substantial improvement in instrumental sta- 
bility since gain and noise fluctuations of amplifiers on 
the two antennas were uncorrelated in time; the only 
correlated signal, that common to the two antennas, was 
from the astronomical source. Communication engineers 
[S E3 had already recognized that correlation tech- 
niques were valuable for retrieving small periodic signals 
in noise, and eventually came to view cross-correlation 
as a method of producing an optimal filter: it selects the 
component of the input signal at one multiplier input 
with waveform equal to the reference signal at the other 
multiplier input |l8j . Viewed in this light, synchronous 
detection is a familiar example of a correlation receiver. 

Adding a four-port circuit (a "hybrid") to correlation 
detection combines signals from two regions of the fo- 
cal plane of a single aperture telescope and redistributes 
them before amplification. This allows the correlation 
technique to be used for single-dish observations. Fig- 
ureUshows the signal flow through a correlation radiome- 
ter; the components to the right of the hybrid are equiv- 
alent to a spatial interferometer's signal path. The ter- 
minology for correlation detection is unfortunately mud- 
dled. Communications engineers use the term correlation 
receiver to describe what radio astronomers usually think 
of as a spatial interferometer or a synchronous detector, 
while single-dish astronomical instruments take the name 
"correlation receiver" [Tol J llT l or the more apt "continu- 
ous comparison receiver" |j, Q ■ More recently, the same 



2 




tional to G{vj etector ) , or v out oc G ((v 2 s 



;)), plus 



Out 



Figure 1: General model for a continuous comparison ra- 
diometer. The central block represents the input hybrid, with 
input ports numbered 1 and 4 and outputs 2 and 3. General 
noise and gain terms n and g and a phase shift ip complete 
the model. 



architecture has been dubbed the "pseudo-correlation re- 
ceiver," apparently based on a detail of the multiplier 
implementation. 

Continuous comparison detection with a single- 
aperture telescope is the complement to the two-element 
spatial interferometer: an interferometer is sensitive to 
the correlated signals from two different regions of an 
aperture plane, while the continuous comparison ra- 
diometer extracts the uncorrelated part of the signals 
from two different regions of a focal plane. The in- 
put hybrid combines voltages from the source and ref- 
erence positions in the focal plane, v s and v r , with differ- 
ent but known phase shifts before amplification. Cross- 
correlating the signals from the two amplifiers with the 
proper phase shift extracts the power difference between 
the two positions, v ou t oc G ((v 2 ) — (v 2 )), where the an- 
gle brackets denote an average over a time long com- 
pared with the reciprocal of the input bandwidth. The 
single-aperture continuous comparison radiometer has 
signal paths that are as similar as possible, so subse- 
quent amplification and processing operate equally on 
the signals from both inputs. An amplifier gain fluctu- 
ation, for instance, has exactly the same effect on the 
signals from both positions in the focal plane. Corre- 
lated terms, including amplifier noise, average away with 
time as 1/VBt, where B is the predetection bandwidth 
and t is the integration time. With a differential mea- 
surement, gain fluctuations have no large noise term to 
amplify, greatly reducing the excess noise across the spec- 
trum. Excess noise is a particular problem for instru- 
ments with wide bandwidths because the intrinsic mea- 
surement noise is proportional to l/\/~B by the radiome- 
ter equation [l^. Fluctuations add noise in individual 
channels and instrumental structure across the spectrum. 

A true differential measurement can greatly improve 
the stability of a radiometer compared with conventional 
total power measurements. Power amplification in a typ- 
ical heterodyne radiometer is about 10 12 , so it is not 
surprising that the most common limit to the stability of 
radiometric measurements is electronic and optical gain 
instability in time. Consider a total power radiometer 
with gain G, input source voltage v s , and system noise 
voltage v n . Its detector output voltage v ou t is propor- 



the small uncorrelated cross term 2(v s v n ) that averages 
away with time as 1/y/Br. Spatial chopping and dif- 
ferencing between source and reference positions on the 
sky largely eliminates the relatively large noise signal on 
average, but even tiny fluctuations in system gain G or 
noise power (tj 2 ) at the chop frequency can easily be much 
larger than the weak signal, A(G(v 2 )) 3> (v 2 ), and can 
dominate the integrated signal. In contrast with Dicke's 
|19| scheme of sequential switching between astronomical 
and reference signals with a total power radiometer, the 
continuous comparison technique's simultaneous treat- 
ment of signal and reference positions provides a true 
differential measurement that greatly reduces the effects 
of time- variable gain fluctuations. Differencing without 
mechanically changing the optical system can also help 
reduce instabilities induced effects from microphonics or 
changing standing-wave structure that affect some types 
of receivers (e.g. local oscillator power modulation from 
a focal plane chopper or nutating secondary) and other 
low-level effects. 



III. CHOICE OF HYBRID 

Figure ^ defines the noise and gain variables for the 
following system analysis. Generalized noise voltages 
nx.Y and voltage gains gx,Y affect the voltages x and 
y at inputs X and Y. The hybrid's output is a phase- 
shifted mixture of its input voltages. Further equiva- 
lent noise voltages ua,b and voltage gains gA,B follow at 
each output. Noise voltages from components following 
the hybrid will be in phase and in quadrature (denoted 
I and Q) with the signal phase, so the noise terms are 
ri = + jnQ/-\/2, where j = \f—l and the total 

noise power is (n 2 ) = (nj) + (rig). Solving for all noise 
and gain components in Fig. ^ is too messy to clearly 
show the circuit's properties. It is clearer to solve two 
basic cases, one with all gain and noise following an ideal 
hybrid and one with all preceding it. Nonideal effects are 
straightforward to include as modifications to these ideal 
cases. 

Most astronomical continuous comparison radiometers 
incorporate a waveguide magic tee 180° hybrid lllij. 
although at least one has used a 90° quasi-optical hy- 
brid and good branch- line 90° h ybr ids now exist at 
frequencies to hundreds of gigahertz [201 |2l| . The corre- 
lator's output has significantly different behavior for the 
180° and 90° hybrids. A lossless hybrid's scattering ma- 
trix relates its output voltages to its input voltages, with 
ports numbered as in Figure ^ as: 
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Here a and j3 are voltage coupling coefficients with a 2 
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(3 2 = 1. Zeros along the diagonal indicate that the ports 
are perfectly matched, and zeros on the cross-diagonal 
indicate no coupling between isolated ports. (A lossless 
hybrid must have these terms equal to zero to satisfy the 
unitary property of a lossless network S matrix 22].) The 
phase angle 9 in Eq. JQ) may vary arbitrarily in theory, 
with 9 = and ir/2 corresponding to realizable devices: 
the fully asymmetrical 180° and fully symmetrical 90° 
hybrids. Both have 180° phase total shifts between the 
outputs, but the shifts are distributed differently relative 
to the input signals. 

Computing the multiplier output v out is easiest when 
the signals are in complex phasor notation with implicit 
time dependence, v(t) = \v\e^. Then the low frequency 



correlator output is v ou t oc (Ke(vAV B )) where va.b are 
voltages at the multiplier input and the asterisk denotes 
the complex conjugate. 

The most useful case has gain and noise following the 
hybrid. Solving for the multiplier's output shows that 
the interesting correlator power-difference signal v ou t oc 
((x 2 ) — (y 2 )) is largest when a term cos((a~ (b — <P + 9) is 
maximum. Here 9 is the hybrid phase defined in Eq. , 
tp is an additional system phase shift shown in Fig.^ and 
(a,b are the gain phases gA9*B = Ge^^ A ~^ B \ Allowing 
for a phase deviation 5, ideally zero, from the maximum 
difference condition, this cosine term is maximum for tp = 
9 + (a — Cb — S. With that substitution the expressions 
become much simpler, and the multiplier output is 



Vout OC 



a/3 ((x 2 ) - (y 2 )) + ((3 2 - a 2 ) (xy) + ^(n A in B i) + (uaqubq) 
--^j=((xnAi) - (yriBi}j + -^={[xn B i) + (ynAi)^ Gcos(S) 
-([uaqubi) - {nAWBQ)) + -j=((xnAQ) + {yriBQ^J 
-^(Wbq) - {yn A c 



for a 180° hybrid (9 = 0), and 



for a 90° hybrid (6 = -n/2). 
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Ideally, the correlator output contains only the power- 
difference term (x 2 ) — (y 2 )- Other uncorrelated cross- 
terms (e.g. (xtia)) average to zero as l/V Br. Not all 
elements within the cross-terms will be completely uncor- 
related in a practical system, with the correlated portions 
producing offsets at the correlator output. Fluctuations 
in system gain G (Eqs.|2]and|3J) scale the offsets and can 
produce error terms that are large compared with the 
signal. An important goal for a continuous comparison 
radiometer is to keep the multiplier output near zero so 
the influence of gain fluctuations will be small. Mini- 
mizing the number and amplitude of correlator offsets is 
important to this end. 

For ground-based radio astronomy the most important 
difference between the circuits with 180° and 90° hybrids 
is likely to be the response to cross-correlation in the in- 



put signals, (xy) . Atmospheric emission in the telescope's 
near field and noise from the telescope's ohmic loss will 
provide correlated voltages between the two inputs. Sup- 
pressing this term requires either good hybrid amplitude 
balance or a flat system phase: at the correlator out- 
put this signal scales with (a 2 — f3 2 ) for the 180° hybrid 
and ~ sin((5) for the 90° hybrid. A factor of ten sup- 
pression implies a hybrid power imbalance no worse than 
0.45/0.55 (0.87 dB) or a maximum phase error term of 
sin(5.7°). Although it is difficult to build hybrids with 
tight amplitude matching across wide bands, wideband 
hybrids can have good phase flatness [2fll2ll |. If the other 
system components have good phase matching then the 
90° hybrid would be the better choice. 

Rejecting correlated input noise can also be useful 
when a common local oscillator signal (LO) is injected 



4 



before the hybrid. Injecting the LO into both ports with 
zero phase shift and equal amplitude will suppress noise 
in the oscillator's wings at the signal frequency. Injecting 
the LO into only one input port will pump both mixers 
but provides no LO noise rejection. 

In any case, the 90° hybrid circuit always rejects corre- 
lated signals introduced after the hybrid better than the 
180° hybrid circuit. These terms, with form (jiaq^bq) 
and {uaitlbi) , are suppressed by the sin((5) factor. Cor- 
relation in these terms can arise from bias fluctuations 
common to both gain paths or noise from the wings of a 
shared local oscillator (the noise and gain model implicit 
in Fig. ^ is generic and can include frequency conversion 
and multiple amplifiers). 

For a nonideal hybrid, coupling between the output 
ports (2 and 3 in Fig. |2J is another important potential 
source of correlator offset. Noise radiated from, or signals 
reflected by, devices following the hybrid can emerge from 
the corresponding nominally isolated port to produce a 
correlated signal. A circulator following the hybrid can 
reduce the offset by an amount equal to the circulator 
isolation at the cost of adding loss before amplification. 
Lack of isolation between hybrid inputs is likely to be less 
of a problem since the signal reflected from the telescope 
or other optics is likely to be small. Further weak corre- 
lated terms will come from noise power emitted from the 
system and reflected back as an input signal (e.g. a frac- 
tion of ha returns as input signal x to produce a nonzero 
(xua))- This term is suppressed if the pathlength for the 
reflected signal is substantially longer than the correla- 
tion length, I ~ c/B, where c is the speed of light. 

Phase switching is a powerful method for removing cor- 
relator offsets and reducing the effects from nonideal col- 
ored noise. In comparison with amplitude modulation 
(Dicke switching) 180° phase switching is very efficient 
because the full signal amplitude is always present at the 
detector. Modulating the phase difference between the 
arms of a continuous comparison radiometer shifts the 
correlated signal output in frequency by an amount equal 
to the modulation frequency. Synchronous detection re- 
covers the correlated signal while rejecting noise fluctua- 
tions at frequencies other than the modulation frequency. 
With phase switching before amplification, ideally just 
following the hybrid, a judicious choice of modulation 
frequency can remove much of the drift and 1// noise as- 
sociated with amplifier noise fluctuations. Since, to high 
order, any residual offsets or offsets from sources outside 
the phase modulation-demodulation boundaries are in- 
dependent of the telescope's pointing, Dicke amplitude 
modulation by chopping on the sky (optically switching 
between two positions on the sky) and a final synchronous 
demodulation will largely remove the remaining offsets. 

Phase switching does not suppress the direct effects 
of amplifier or multiplier gain or phase fluctuations, but 
symmetry can reduce their influence. Equations (0) and 
show that the product of the voltage gains scales the 
input power difference as G((x 2 ) — (y 2 ))- A small gain 
fluctuation common to both chains introduces an am- 



plitude error to the difference signal at the correlator 
output, but adds no error signal when the offsets are 
negligible. Multiplier gain fluctuations have the same 
effect as amplifier gain fluctuations in this case. The 
situation is different when the amplifier gain fluctua- 
tions are differential-mode instead of common-mode for 
the two amplifier chains. Faris jjj calculates the ef- 
fect of a varying gain imbalance between the two chains 
for gz{f,t) = [1 + a(t)]5i(/), where 5i, 2 (/) are the 
complex amplifier voltage gains of the two chains and 
a(t) is a zero-mean random variable that describes the 
differential-mode fluctuations. Fluctuations increase the 
output variance by a factor of (1 + (a 2 (t))) compared 
with the case of purely common-mode gain fluctuations 
between the two arms. Similar effects arise from differen- 
tial phase fluctuations between the two arms. Strategies 
for minimizing differential-mode fluctuations could in- 
clude biasing the amplifiers from a common power supply 
and keeping good thermal contact between corresponding 
parts of each chain. 

The second case, with amplification preceding the hy- 
brid's loss, is the obvious choice for maximizing the re- 
ceiver sensitivity but negates much of the continuous 
comparison architecture's advantage. The maximum dif- 
ference signal in this case is 
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for both 180° and 90° hybrids. For this configuration to 
be useful nearly exact matching of both the noise power 
and power gains would be necessary. Slight imbalances in 
loss and gain will produce large offsets at the correlator 
output. 



IV. SPECTROSCOPY WITH AN ANALOG LAG 
CROSS-CORRELATOR 

Sensitive spectroscopy (multichannel radiometry) over 
broad bandwidths is important for observations of wide 
spectral lines from distant galaxies, searches for lines at 
with uncertain frequency, and measurements of pressure- 
broadened lines in planetary atmospheres. Spectrometer 
bandwidths can be tens of gigahertz with channel band- 
widths of tens of megahertz. Such broad bandwidths 
place stringent requirements on system stability, so it is 
natural to pair wideband spectrometers with the contin- 
uous comparison architecture. Spectroscopy with a cor- 
relation radiometer requires measurement of the cross- 
correlation function over a range of time lag, or delay. 
Analog lag correlators use purely analog components to 
obtain the cross-correlation function Rab{t) as a func- 
tion of lag r: 



1 f T 

Rab (r) = hm — / v A (t) -v B (t + T)dt 

I -^OCl Zl J _ rp 



(5) 
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Tapped transmission lines provide the time delays r, 
transistor multipliers form the product of the two input 
voltages VAif) and Us(t+r), and low-frequency electron- 
ics integrate the multiplier output to provide the time 
average. 

A Fourier transform of the cross-correlation function 
yields the power density cross-spectrum. Transform- 
ing the correlation function to recover the spectrum is 
slightly more complicated than making a direct Fourier 
transform for analog correlators because the signal is not 
sampled at perfectly regular intervals. Although the me- 
chanical spacing between the microwave signal taps along 
the transmission line is well defined by the traces on the 
circuit board, frequency-dependent component variations 
cause small erratic variations in the electrical delays be- 
tween multipliers. A direct Fourier inversion to find the 
spectrum is not possible because the transform kernel's 
phase term cannot be reduced to a separable product of 
the delay time and frequency, as might be the case for 
simple line dispersion It is possible to correct the 
irregular sampling in software by measuring the spec- 
trometer's response to a series of monochromatic signals 
at known frequencies, then expanding the astronomical 
input signal on these measured functions 0. 

An interesting effect of this calibration scheme is that 
it defines the spectrometer's internal phase: by defini- 
tion signals in phase with the calibration are real, and 
those with a relative 90° phase shift are imaginary. This 
property can be used to eliminate one of the phases in 
the usual complex cross-correlation measurement. When 
the calibration signals are injected at the radiometer's in- 
put the calibration and measured signals share the same 
phase shifts through the entire instrument, so the mea- 
sured signal is purely real; the imaginary component 
contains only noise. A single cross-correlator can there- 
fore measure the cross-correlation function. While purely 
real cross-correlations are unusual in most spectrometers, 
there is no fundamental reason that they cannot exist. A 
real correlation function has the convenient property of 
even symmetry in the lag domain, so the positive (or neg- 
ative) lags alone contain all the necessary information to 
recover the spectrum. Although it is possible to build 
a full complex correlator and calibrate it at lower fre- 
quency, injecting the phase calibration signals at the in- 
put to the entire radiometer yields spectra at full spectral 
resolution with half the number of lags. This is not only 
a significant savings in spectrometer cost and complexity, 
but eliminates requirements on phase-matching between 
the two amplifier and processing chains. Digital corre- 
lators do not readily share this property because their 
topologies give them an intrinsic symmetry and phase 
related to the position of the zero-lag multiplier. A full 
complex correlator, or close phase matching across the 
receiver band, is needed for spectroscopy with a digital 
system and direct transform. 

A simple laboratory continuous comparison radiome- 
ter, sketched in Figure [3 verifies that a WASP2 (Wide- 
band Astronomical SPectrometer) analog lag correlator 




Figure 2: Block diagram of laboratory continuous comparison 
receiver test setup. 
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Figure 3: Spectra from the two inputs of a laboratory cor- 
relation radiometer with a WASP2 spectrometer configured 
as a cross-correlator. The dotted line is a network analyzer 
measurement of the filter transmission. 



I] properly produces power difference cross-spectra with 
this calibration scheme. The hybrid for the experiment 
was an off-the-shelf stripline 90° 2-4 GHz device. Cables 
lengths between the hybrid and cross-correlator brought 
the zero time-lag position close to one end of the multi- 
plier ladder, for maximum spectral resolution, but were 
not otherwise trimmed for length or phase matching. A 
broadband noise diode and 300 MHz filter generated an 
artificial spectral line that could be connected to either 
hybrid input, denoted X and Y in Fig. [21 Phase calibra- 
tion signals were fed into input X alone. Figure [3] shows 
that the spectrometer works as predicted. The artificial 
line at input X, with Y terminated, produces a positive 
spectral line, and the same signal at input Y, with X ter- 
minated, produces a negative spectral line. After remov- 
ing the noise source's intrinsic spectral shape by divid- 
ing the raw spectra by spectra of the noise source alone, 
the filter center frequency, shape, and loss matched net- 
work analyzer measurements (Fig.|3J). This confirms that 
the cross-correlation function is purely real. The sum of 
the two spectra is very close to zero. Residual phase er- 
rors scatter power across the spectrum at about 0.5% of 
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the peak line intensity, a well-understood dynamic range 
limit rather than a noise offset pj. This fixed-pattern 
structure subtracts well with beamswitching. 

These experimental results show that the combination 
of the continuous comparison architecture with a ana- 
log lag correlator is a very promising method for spec- 
troscopy over very wide bandwidths. 
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